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ON QUANTUM COMPUTATION, ANYONS, AND 

CATEGORIES 


ANDREAS BLASS AND YURI GUREVICH 

Abstract. We explain the use of category theory in describing 
certain sorts of anyons. Yoneda’s lemma leads to a simplification 
of that description. For the particular case of Fibonacci anyons, 
we also exhibit some calculations that seem to be known to the 
experts but not explicit in the literature. 


1. Introduction 

This paper attempts to explain the use of category theory in de¬ 
scribing certain sorts of anyons. These are rather mysterious physical 
phenomena which, one hopes, will provide a basis for quantum com¬ 
puting needing far less error correction than other approaches. 

The first author of this paper has long been a fan of category theory; 
even as a graduate student, he was described by one of his professors 
as “functorized”. The second author has been far more skeptical about 
the value of category theory in computer science, because of its distance 
from applications and because of the peril of potential (and in some 
cases actual) over-abstraction. In 2012, both authors began working 
with the Quantum Architectures and Computing (QuArC) Group at 
Microsoft Research and found anyons to be near the top of the group’s 
agenda. Seeing calculations and applications that use unitary matrices 
to represent braiding of anyons, we naturally wondered what Hilbert 
space these matrices are intended to operate on. We made rather a 
nuisance of ourselves by asking different people, on different occasions, 
what anyons actually are, from a mathematical point of view. Are 
they Hilbert spaces? Are they vectors in a Hilbert space? Are they 
something else? It turned out that the only mathematically sound 
answer in the literature involved a special sort of categories, modular 
tensor categories^ So the second author agreed that categories can be 
quite relevant to important applications in computer science. 

^Other answers explained the physics, in terms of excitations, but these matters 
are not the subject of this paper, which is specifically about mathematics except 
for the introductory material summarized in Section [2] 
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Our purpose in this paper is to describe some of the ideas surround¬ 
ing categories and anyons in general and the special case of Fibonacci 
anyons and their category description. We hope that our presentation 
will be accessible and useful for mathematicians and computer scien¬ 
tists who have some acquaintance with the basics of category theory. 
Where we need to go beyond the basics, we explain, albeit briefly, the 
concepts from category theory that we use. We have also included a 
section describing the physical background that this mathematics is 
intended to formalize. 

To describe more of our motivation for studying anyons, we need 
to presuppose some general information that will be explained in later 
sections of this paper. In particular, we shall refer to the fusion rule 
r 0 r = r © 1 for Fibonacci anyons r (and the vacuum 1). We hope 
that the following paragraphs will give the reader a rough idea of what 
we are looking at, and that re-reading them after the rest of the paper 
will provide a less rough idea. 

In contrast to what occurs elsewhere in quantum theory, the states 
(represented, as usual, by vectors in Hilbert spaces, up to scalar mul¬ 
tiples) in the modular tensor category picture are ways in which one 
conhguration can fuse to form another configuration]^ They are not the 
configurations themselves. For example, in the Fibonacci case, there is 
a 2-dimensional Hilbert space of ways for three anyons to be regarded 
as (or to fuse into) one anyon; this is the Hilbert space Hom(r(8)r(8)r, r). 

When we first heard about Fibonacci anyons, we thought that the 
fusion rule r © r = r © 1 meant that, if we put two r anyons together, 
then the result might look like one r anyon or like the vacuum (this 
much is true in the modular tensor category model) and that the gen¬ 
eral result would be a superposition of these two alternatives. But 
the model doesn’t allow such superpositions. Nor does the model say 
anything about the probabilities of the two possible outcomes. 

Instead, we get superpositions of the following sort. Start with three 
r’s. Fuse the first two to get one r or vacuum. If you got vacuum, then 
the overall result is one r, namely the third of the original ones, which 
you haven’t yet fused. If, on the other hand, fusing the first two r’s 
gives a r, then fusing that with the third r might produce a r. (It might 
also produce vacuum, but that’s irrelevant for the present discussion.) 
So we have two ways to end up with one r, according to whether the 
first two r’s fused to vacuum or to r. And it is these two ways that 
the model allows superpositions of. Another possibility for getting two 
ways here is to fuse the last two r’s first and then fuse the result with 

^For more on the notion of fusion, see Remark [1] at the end of this introduction. 
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the first r. These two form another basis of the same 2-dimensional 
Hilbert space of “ways”. The relation between the two ways is (part 
of) the associativity isomorphism of the modular tensor category. Yet 
another possibility would begin by fusing the hrst and third r’s. The 
modular tensor category representation of this possibility would use a 
braiding isomorphism to move the first anyon to be adjacent to the 
third (or vice versa), and it would depend on the path along which 
that anyon is moved around the second one. 

In Section |21 we give a general introduction to anyons from the point 
of view of physics and quantum computation. That section is intended 
to give the reader a rough idea of what anyons are and why researchers 
in quantum computation would be interested in them. The treatment 
here is quite superficial, and we give references for more detailed treat¬ 
ments. 

In Section [3l we gradually introduce modular tensor categories, and 
we explain how they are intended to be used to describe anyons. This 
section borrows heavily from the axiomatization given in [5], but with 
some modihcations and rearrangements. 

Section 0] is devoted to an application of one of the central theo¬ 
rems of category theory, known as Yoneda’s Lemma, to producing a 
simplihed view of modular tensor categories. 

Finally, in Section |5l we consider the special case of Fibonacci 
anyons. This special case is unusually simple in some respects. Never¬ 
theless (or perhaps therefore) it occupies a prominent place in quantum 
computing research. Section [5] begins with a general description of Fi¬ 
bonacci anyons and then exhibits some calculations, whose results seem 
to be well known to some in the quantum computing community but 
which we have not been able to find written down in the literature. 

More detailed treatments of modular tensor categories are available 
in the papers [9] of Panangaden and Paquette and m of Wang. Much 
of our exposition is based on the former. For other aspects of anyons 
and topological quantum computation, see, for example, [5] and the 
references there. 

Remark 1. We encountered numerous explanations of the notion of 
fusion of anyons, and they seemed to contradict each other. At one 
extreme was the picture of fusion as a physical process in which anyons 
are brought into spatial proximity with each other and energy is re¬ 
leased as they form a new anyon (or perhaps annihilate each other). A 
minor modihcation of this picture is that energy need not be released; 
it might actually be consumed in the process. Another picture, how¬ 
ever, did not insist that the anyons be brought together. They could 
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remain far apart, and a suitable global measurement of the system’s 
quantum numbers could reveal how they “fused”. A path to reconciling 
these apparently contradictory pictures is suggested by a comment at 
the end of Section II.A of [8]; the idea is as follows. Consider several 
anyons, which we intend to fuse. As long as they are far apart, the 
various possible results of their fusion have energies that are very close 
together. (In technical terms, the ground state of the system is very 
nearly degenerate.) So the different fusion results can be distinguished 
in principle but not practically. When the anyons are brought closer 
together, though, the energy differences between the fusion possibilities 
become larger, and so it becomes practical to distinguish these possi¬ 
bilities. Thus, the discrepancy between various views of fusion seems 
to be largely a discrepancy between what can be observed in principle 
(or what is “really” happening) and what can be detected in practice. 

2. Quantum theory and anyons 

This section is a superhcial summary of a small part of quantum the¬ 
ory and some basic information about anyons. The physics described 
here is intended merely to provide an orientation for understanding the 
mathematics in the rest of the paper. 

2.1. Quantum Mechanics. In quantum theory, the state of a phys¬ 
ical system is typically represented by a non-zero vector in a complex 
Hilbert space "H, but all non-zero scalar multiples of a vector represent 
the same state. Thus, the states constitute the projective space asso¬ 
ciated to 'H. Because of the freedom to adjust scalar factors, one often 
imposes the normalization that the vectors representing a state should 
have norm 1; there still remains a freedom to adjust the phase, i.e., a 
scalar factor of absolute value 1. 

If a system has an observable property with inhnitely many possible 
values, for example position or momentum, then the Hilbert space of 
its states must be inhnite-dimensional. In quantum computing, how¬ 
ever, one usually ignores many such properties and concentrates on 
only a small number (often only one) of properties with only hnitely 
many possible values. As a result, one deals with hnite-dimensional 
Hilbert spaces. (This simplihcation is analogous to modeling a classi¬ 
cal computer by a conhguration of bits, not taking account of its other 
physical properties, like position or momentum or temperature, unless 
these threaten to interfere with the bits of interest.) 

The automorphisms of a Hilbert space "H are the unitary trans¬ 
formations, i.e., the linear bijections that preserve the inner product 
structure. These play several important roles, both in physics and in 
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quantum computation. First, they provide the dynamics of isolated 
quantum systems. That is, the state of an isolated system will evolve 
in time by the action of a one-parameter group (the parameter being 
time) of unitary operators^ Second, if a system has symmetries, i.e., 
if it is invariant under some transformations, then these transforma¬ 
tions are usually modeled by unitary operators]^ Finally, the design 
of quantum algorithms is based on unitary operators. We want the 
system to evolve from a state that we know how to produce to another 
state from which we can extract useful information by a measurement. 
That evolution is described by a unitary operator. So an algorithm de¬ 
signer wants to hnd unitary operators that represent a useful evolution 
of a state. In addition to hnding such operators, we want to represent 
them as compositions of simpler ones, called gates, that we know how 
to implement. 

Where classical computation uses bits, whose possible values are de¬ 
noted by 0 and 1, quantum computation uses qubits. A measurement 
of a qubit produces two possible values; the qubit itself is represented 
by a 2-dimensional Hilbert space, in which a certain orthonormal basis, 
usually written {|0), |1)}, corresponds to the two values. In contrast to 
the classical case, though, the Hilbert space structure provides many 
other states in addition to these two basic ones. Any non-zero linear 
combination of |0) and |1) represents a possible state of the system. If 
the state is represented by the unit vector a:|0) -|- i/|l), then measuring 
the qubit in the {|0), |1)} basis will produce the outcome 0 with prob¬ 
ability |xp and the outcome 1 with probability ||/p. Such a state is a 
superposition of the two basic states. More precisely, this state vector 
is the superposition, with coefficients x and y, of the vectors |0) and 
|1), respectively. 

It is more accurate to speak of superposition of vectors than of su¬ 
perposition of states. The reason is that, although phase factors don’t 
affect the state represented by a vector, relative phases do affect super¬ 
positions. Thus, for example, although |1) and —11) represent the same 
state of a qubit, the superpositions (|0) -|- |l))/\/2 and (|0) — |l))/\/2 
represent quite different states. 

It is almost true in general that, for any two states of any quantum 
system, any superposition of the associated vectors also represents a 
possible state of that system. The word “almost” in the preceding 

^Here we use the so-called Schrodinger picture of quantum mechanics. A physi¬ 
cally equivalent alternative view, the Heisenberg picture, has the states remaining 
constant in time, while the operators modeling properties of the state evolve by 
conjugation with a one-parameter group of unitary operators. 

few discrete symmetries can be modeled by anti-unitary transformations. 
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sentence refers to the possibility of supers election rules. These rules 
specify that, for certain quantities, like electric charge, it is impossible 
to superpose two states with different values of those quantities. Thus, 
when discussing a system for which several values of the electric charge 
can occur, we are, in effect, dealing with several separate Hilbert spaces, 
called superselection sectors, one for each value of the charge. One can, 
and sometimes one does, form the direct sum of these Hilbert spaces to 
obtain a Hilbert space containing all the possible states of that system, 
but most of the vectors in that direct sum, involving superpositions 
with different charges, do not represent physically possible states. We 
prefer, in this paper, to deal with superselection sectors as separate 
Hilbert spaces and forgo their direct sum. For more information about 
superselection rules, see [1]. 

In reality, there are very few superselection rules — arising from 
certain conserved quantities like electric charge, baryon number, and 
parity — but in the study and application of anyons one often artifi¬ 
cially adds superselection rules, and we shall encounter such rules in 
the category-theoretic treatment below. This amounts to deciding not 
to consider superpositions of vectors from certain Hilbert spaces, i.e., to 
consider those superselection sectors separately rather than considering 
their direct sum. 

In the presence of superselection rules, the operators that one consid¬ 
ers are operators acting on each of the superselection sectors separately. 
In the case of true superselection rules, the dynamics of the system and 
any gates that one could construct are given by unitary operators act¬ 
ing on each sector separately. In the case of artificial superselection 
rules, nature may not cooperate with our artihcial rules, and states in 
one sector may evolve out of that sector. Such evolution interferes with 
our understanding and intentions; it is often called “leakage” and one 
strives to avoid it. 

In addition to the unitary operators mentioned above, Hermitian 
(or self-adjoint) operators on the Hilbert space of states also play an 
important role in quantum mechanics, because they model observable 
properties of a system. The connection between Hermitian operators 
and (real-valued) observables is easy to describe in the case of hnite- 
dimensional Hilbert spaces hE Let the Hermitian operator A have 
(distinct) eigenvalues Oi,..., a^, with associated eigenspaces Si,, Sk. 
(Some of these eigenvalues may have multiplicity greater than 1, but 


®In the infinite-dimensional case, the description is similar but one must take 
into account the possibility of a continuous spectrum of the operator, in addition 
to or instead of discrete eigenvalues. 
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they are to be listed only once among the ads. The associated Si will 
then have dimension greater than 1.) These eigenspaces are orthogonal 
to each other, and their sum is all of "H. Any unit vector I'ljj) G Ti can 
be expressed as the sum of its projections |<pj) to the subspaces S'*. 
Measuring A on a system in state produces one of the eigenvalues 
tti] the probability of getting the result a* is the squared norm of the 
projection, |||v7i)ll^- Note that the dimension of "H is an upper bound 
for the number of distinct eigenvalues a* of any Hermitian operator on 
'H. In particular, any measurement performed on a qubit will have 
at most two possible outcomes. It is in this sense that a qubit is the 
quantum analog of a classical bit. 

2.2. Anyons. To understand anyons, it is useful to recall first that or¬ 
dinary particles are of two sorts, bosons and fermions. These differ in 
several respects, beginning with the action of spatial rotations on the 
corresponding Hilbert spaces. For particles in ordinary 3-dimensional 
space, the group SO(3) of Euclidean rotations of that space acts on 
the states of the particle. (More precisely, the group of all Euclidean 
motions acts, but we abstract from the particle’s position and consider 
only its orientation in space; thus we ignore translations and consider 
only the group of rotations.) Because the vector representing a state 
is dehned only up to a phase factor, the action of the rotation group 
is not a representation in the usual sense but a projective representa¬ 
tion. This means that each rotation g of physical 3-dimensional space 
is represented by a unitary operator p{g) on the Hilbert space, but this 
p{g) is unique only up to a phase factor. It is customary to make some 
arbitrary choice of these phase factors, so that we can speak unam¬ 
biguously of p{g). The arbitrariness of the choice is, however, reflected 
in the fact that p{gh) and p{g)p{h) need not be equal but can differ 
by a phase factor. Furthermore, p and p' are considered equivalent 
representations if they differ only by these arbitrary phase factors. It 
is reasonable to ask, in this connection, why the operators p{g) need to 
be unitary or even linear, rather than only linear up to phase factors. 
The reason is that, unlike absolute phases, relative phases are relevant 
in superpositions, so physical symmetries must preserve them. 

It turns out that any projective representation p of SO{3) is given 
by a genuine unitary representation p of the universal covering group 
of ^(^(d), namely SU{2) (see for example [1] and [lO]). That is, if 
p : SU{2) —)■ SO{3) denotes the 2-to-l projection map, we have pop 
equivalent to p. More concretely, it means that there are two sorts 
of projective representations of SO{3), up to equivalence. One sort is 
the ordinary unitary representations of S'0(3); the other is given by 
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unitary representations of SU (2) that send the non-trivial element —I 
of the kernel of p to the operator —I. (Throughout this paper, we 
use I, sometimes with subscripts, to denote identity transformations, 
functions, morphisms, etc.) The first sort of representation corresponds 
to bosons, whose state vectors (not merely their states) are unchanged 
when rotated gradually through a full revolution. The second sort 
corresponds to fermions, where a rotation through 27r changes the state 
vector by a sign. 

A second distinction between bosons and fermions, even more im¬ 
portant for our purposes, is the behavior of systems of several iden¬ 
tical particles. Because the particles are identical, any permutation 
of the particles leaves the state unchanged and therefore changes the 
state vector by at most a phase factor. As a result, we have a one¬ 
dimensional projective representation of the symmetric group. Again, 
it turns out that there are just two possibilities (both of which are 
actual unitary representations of the symmetric group). Either all per¬ 
mutations leave the state vectors unchanged, or the even permutations 
leave the state vectors unchanged while the odd permutations reverse 
the vectors’ signs. 

A deep theorem of relativistic quantum field theory, the spin- 
statistics theorem, says that these two behaviors of multi-particle states 
under permutations exactly match the two behaviors of single-particle 
states under rotations. Interchanging two identical bosons leaves the 
state vector of the pair unchanged; interchanging two identical fermions 
reverses the sign of the state vector. 

The preceding discussion of bosons and fermions depends crucially on 
the fact that the particles are in ordinary 3-dimensional space. If par¬ 
ticles were conhned to a 2-dimensional space, more possibilities would 
arise. 

Specihcally, the rotation group in two dimensions, SO{2), has more 
sorts of projective representations than SO(3) does; the reason is ulti¬ 
mately that the universal covering group of the circle group SO{2) is 
the additive group of real numbers, and the covering projection is not 
2-to-l but cx)-to-l. The result is that a gradual rotation of a particle 
through 27r can multiply its state vector by an arbitrary phase factor, 
not just ±1. The possibility of getting any phase here led to the name 
anyon. 

Reducing the dimensionality of space from 3 to 2 also affects the 
possibilities for permuting identical particles. For simplicity, consider 
the case where there are just two particles, and we interchange them. 
We can perform the intechange gradually, in the plane, by rotating the 
2-particle system counterclockwise by vr around the midpoint between 
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the particles. Alternatively, we can achieve the same interchange by 
a clockwise rotation. In 3-dimensional space, these two options are 
equivalent in the sense that they can be gradually deformed into each 
other, by rotating the plane of the particles’ motion about the line 
through the particles’ initial positions. In 2-dimensional space, there 
is no such deformation without making the particles collide. Winding 
one particle around the other any number of times, we get inhnitely 
many ways to achieve one and the same permutation. With more 
than two particles, there are more complicated ways to achieve the 
same permutation by moving the particles around in the plane. As 
a result, in place of (projective) representations of symmetric groups, 
we have representations of braid groups. For example, in the case of 
two particles, in place of the group of two possible permutations of the 
particles, we have the group of all integers, with integer n representing a 
counterclockwise rotation by nvr (and negative n representing clockwise 
rotations). 

The preceding discussion was oversimplified in that (among other 
things) when moving particles around each other, we ignored any ro¬ 
tation that the individual particles might have undergone during the 
motion. A more accurate presentation would need to suitably combine 
the braid and rotation groups. 

2.3. Anyons in Reality. As explained above, anyons do not occur in 
3-dimensional space; it is necessary to reduce the number of spatial di¬ 
mensions to 2. Since we live in a 3-dimensional space, will we ever find 
anyons? It turns out that anyon-like behavior occurs for certain excita¬ 
tions in materials that are so thin as to be effectively two-dimensional. 
A detailed discussion of this would take us too far from the purpose of 
this paper, so we refer the reader to Section 1.1 of [S]. 

We emphasize, however, that the anyons are not what one would or¬ 
dinarily think of as “particles” but rather excitations in some medium, 
which exhibit particle-like behavior. It should be noted in this connec¬ 
tion that it is not unusual, in other contexts, for excitations to behave 
like particles and thus to be analyzed mathematically as if they were 
particles. For example, vibrational excitations in crystal lattices are 
treated as particles called phonons. Similarly, photons are excitations 
of the electromagnetic field. In quantum field theory, all particles are 
excitations of the corresponding fields. 

2.4. Anyons in Quantum Computation. Quantum computation 
is unpleasantly susceptible to environmental disturbances. Its advan¬ 
tages over classical computation depend on maintaining superpositions 
of state vectors, with high precision in the coefficients of those vectors. 
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Small disturbances can easily modify those coefficients or, indeed, de¬ 
stroy superpositions altogether. Signihcant effort must therefore be de¬ 
voted to error correction, and this makes algorithms slower and harder 
to design. 

It has been suggested [6] that qubits could be more robust, i.e., less 
susceptible to disturbances, if they were implemented using certain 
sorts of anyons. For example, if qubits were encoded in the way two 
anyons wind around each other, then this winding, being a topological 
property of the system, would be robust. A small disturbance in the 
actual motion of the anyons would leave the winding number intact. 
This hope of reducing the error correction needs of quantum computing 
has motivated much of the current interest in anyons. 

In this approach to quantum computation, braiding of anyons serves 
not only to store information but also to process it. In general, as men¬ 
tioned above, quantum computation proceeds by initializing a quan¬ 
tum state, then applying a unitary transformation to it, and hnally 
measuring some observable in the resulting transformed state. The 
unitary transformation used here must be designed so that a feasible 
measurement produces a useful result. In addition, there must be a 
way to implement the unitary transformation as the composition of a 
sequence of simpler unitary tranformations, usually called gates in this 
context. In the anyon approach to quantum computation, the most ba¬ 
sic unitary gates arise from the braiding of anyons around each other, 
and a crucial question is whether these gates are universal in the sense 
that arbitrary gates can be approximated by composing the basic ones. 

It is worth noting explicitly that, in this picture, a qubit is not 
encoded in the state of a single anyon but rather in a whole system of 
several anyons. This feature will be quite prominent in the category 
picture described in the rest of this paper. 

3. Modular tensor categories 

In this section we describe the category-theoretic structure that has 
been developed to support a mathematical theory of anyons. Much of 
what we describe here is in [S], though we have modihed some aspects 
and rearranged others. 

Throughout this section, we let M be a category, intended to describe 
the quantum-mechanical behavior of a system of anyons. A will carry 
several sorts of additional structure, roughly classihed as “additive” 
and “multiplicative” structure, all subject to various axioms. We de¬ 
scribe the structures and the axioms a little at a time. We begin with 
the additive structure, because this is where Hilbert spaces enter the 
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picture, so it is the basis for the connection with the usual formalism 
of quantum theory. 

The vectors in our Hilbert spaces will be the morphisms of A. Specif¬ 
ically, for each pair of objects X, Y of A, the set Hom(X, Y) of mor¬ 
phisms from X to Y will have the structure of a Hilbert space. So 
we have many Hilbert spaces, one for each pair X, Y of objects. Some 
of these Hilbert spaces will be mere combinations of others, but there 
will still be several different “basic” Hilbert spaces. This means physi¬ 
cally that we regard the system as being subject to superselection rules, 
which keep these Hilbert spaces separate. 

We assume familiarity with some basic notions of category theory, 
specifically, the notions of product (including terminal object, which 
is the product of the empty family), coproduct (including initial ob¬ 
ject), equalizer, coequalizer, monomorphism, epimorphism, isomor¬ 
phism, functor, and natural transformation. Definitions and examples 
can be found in or [3l Chapter 1]. 

3.1. Additive Structure. We begin by requiring A to be an abelian 
category. This requirement, formulated in detail below, provides a well- 
behaved addition operation on each of the sets Hom(X, H), although 
the requirement is formulated in purely category-theoretic terms and 
does not explicitly mention this addition operation. 

Axiom 1 (Abelian). A is an abelian category. That is 

(1) There is an object 0 that is both initial and terminal. A mor¬ 
phism that factors through this zero object will be called a zero 
morphism and denoted by 0. Note that each Hom(X, Y) con¬ 
tains a unique zero morphism. 

(2) Every two objects have a product and a coproduct. 

(3) For every morphism a : X ^ Y, the pair a, 0 has an equalizer 
and a coequalizer. These are called the kernel and cokernel of a. 

(4) Every monomorphism is the kernel of some morphism, and 
every epimorphism is the cokernel of some morphism. 

This axiom has a surprisingly rich collection of consequences, de¬ 
veloped in detail in Chapter 2 of [3]. We list here only some of the 
highlights, which will be important for this paper, and we refer the 
reader to [3] for the proofs and additional information. 

Proposition 1 ([3], Theorem 2.12). Any morphism that is both monic 
and epic is an isomorphism. 

(More generally, as one can easily check, in any category, any equal¬ 
izer that is an epimorphism is an isomorphism.) 
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Proposition 2 (j3], Theorem 2.35). The product and coproduct of any 
two objects coincide. 

That is, given two objects X and T, there is an object X (BY that 
serves simultaneously as the product of X and Y, with projections 
px ■ X (BY X and py ■ X (BY —^Y, and as the coproduct of X 
and y, with injections ux ■ X ^ X (B Y and uy ■ Y —>■ X (B Y. (If 
X = Y^ then our notations for the projections and injections become 
ambiguous, and we use Pi,P 2 ,Wi,M 2 instead.) For brevity, we often 
refer to X © T as the sum of X and P, rather than as the product or 
coproduct. 

As a product, X(BX admits a diagonal morphism Ax : X ^ X©A, 
namely the unique morphism whose composites with both projections 
are the identity morphism Ix of X. Dually, as a coproduct, it admits 
the folding morphism Vx '■ X (B X ^ X, whose composites with both 
of the injections are Ix- Using the diagonal and folding morphisms, we 
can define a binary operation, called addition, on Hom(X, Y) for any 
objects X and Y. Given f,g:X^Y,we define f + g : X ^ Y to be 
the composite 



where / © is obtained from the functoriality of products (or of co¬ 
products — they yield the same result). 

Proposition 3 ([3], Theorems 2.37 and 2.39). This addition opera¬ 
tion makes each Hom(X, Y) an abelian group, with the zero morphism 
serving as the identity of the group. Composition of morphisms is ad¬ 
ditive with respect to both factors; that is, when either factor is fixed, 
the composite f o g is an additive function of the other factor. 

Axiom 2 (Vectors). Each of these abelian groups Hom(X, V) car¬ 
ries an operation of multiplication by complex numbers, making 
Hom(X, Y) a vector space over C, and making composition of mor¬ 
phisms bilinear over C. 

The complex vector spaces Hom(X, Y) will play the role of quantum- 
mechanical state spaces. For this purpose, they should also be equipped 
with inner products, making them Hilbert spaces, but, following [9], we 
refrain from assuming an inner product structure at this stage of the 
development]^ It turns out that much of what we shall do later does 


®In fact, inner products are never explicitly assumed in [^. They are, however, 
implicit in the statement, in Section 5.1 of [5], that certain bases “are - of course 
- related by a unitary trairsformation”. 
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not depend on the availability of inner products in the vector spaces 
Honi(X,r). 

An object S in the abelian category A is called simple if S' ^ 0 
and every monomorphism into S is either a zero morphism or an iso¬ 
morphism. In other words, S' is a non-zero object with no non-trivial 
subobjects. Because of the abelian structure of A, this dehnition can 
be shown (using [3l Theorem 2.11]) to be equivalent to its dual: A 
non-zero object is simple if and only if it has no non-trivial quotients, 
i.e., every epimorphism out of S' is either a zero morphism or an iso¬ 
morphism. 

Axiom 3 (Semisimple). Every object in A is a hnite sum of simple 
objects. 


This axiom considerably simplihes the structure of the vector spaces 
Hom(A, y). In the hrst place, as shown in [SJ Section 2.3], morphisms 
from a sum 0^. Sj to another sum S(, are given by matrices of 
morphisms between the summands. Specihcally, the matrix associated 
S'(, has as its a, b entry the composite 


Si ^ 0s, ^ 0s; 

j k 


SI. 


Composition of morphisms in A corresponds to the usual multiplication 
of matrices. 

Furthermore, when the summands are simple, we have the following 
additional information about the matrix entries, a generalization of 
Schur’s Lemma in group representation theory. 


Proposition 4. If f : S ^ S' is a morphism between two simple 
objeets, then f is either the zero morphism or an isomorphism. 


Proof. The kernel of / is a monomorphism into S', and if it is an iso¬ 
morphism then / is zero. So, by simplicity of S', we may assume that 
the kernel of / is zero and therefore (by [5|, Theorem 2.17*]) / is a 
monomorphism. Similarly, by considering the cokernel of / and invok¬ 
ing the simplicity of S', we may assume that / is an epimorphism. But 
then, by Proposition [H / is an isomorphism. □ 

The last axiom in this subsection combines two hniteness assump¬ 
tions. 


Axiom 4 (Finiteness). (1) There are only hnitely many non¬ 
isomorphic simple objects. 

(2) Each of the vector spaces Hom(X, Y) is hnite-dimensional over 


C. 
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The first of these two hniteness requirements is merely a technical 
convenience. The second, however, gives the following important in¬ 
formation about the endomorphisms of simple objects. 

Proposition 5. If S is a simple object, then Hom(S', S') = C. 

Proof. The operation of composition of morphisms is a multiplication 
operation that makes the vector space Hom(S', S) into an algebra over 
C. Since S is simple. Proposition 0] says that every non-zero element of 
this algebra is invertible. That is, Hom(S', S) is a division algebra over 
C. But C is algebraically closed, so the only hnite-dimensional division 
algebra over it is C itself. □ 

Note that the isomorphism Hom(S', S') = C in this proposition can 
be taken, as the proof shows, to be an isomorphism of algebras, not just 
of vector spaces. In particular, the identity morphism of S' corresponds 
to the number 1. 

Combining this proposition with our earlier observations about ma¬ 
trices, we hnd that any morphism / ; ®, Sj ^ S'^ between any two 

objects in A is given by a matrix whose entries are complex numbers. 
Moreover, the a, b entry is 0 unless St, = S'^. From this observation, 
it easily follows that, when an object X of M is expressed as a sum 
0^ Sj of simple objects, the isomorphism types of the summands Sj 
and their multiplicities are completely determined by X. That is, the 
representation of X as a sum of simple objects is essentially unique. 

3.2. Multiplicative Structure. In this subsection, we introduce the 
multiplicative structure that makes A a braided monoidal category. 
The central idea is that, if objects X and Y represent certain anyons, 
then X®Y should represent a system consisting of both of these anyons. 
We must, however, remember that the Hilbert spaces that occur in this 
context are not the objects of A but the vector spaces of morphisms 
between the objects. 

A system consisting of two anyons of types X and Y would, if mea¬ 
sured as a whole, appear as another anyon, whose type might not be 
entirely determined by the types X and Y. Formally, this means that 
X ®Y is a sum of several simple objects. Furthermore, there might be 
several “ways” for a composite system to appear as having a particu¬ 
lar type Z, modeled as several morphisms from X ®Y to Z, and our 
Hilbert spaces will also contain superpositions of these. 

The multiplicative structure will also include a unit object 1; its 
intended interpretation is the vacuum. Thus, 1 ® X and X (8) 1 amount 
to just X because a system consisting of X and nothing is the same as 
X. 
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The first aspect of multiplicative structure can be stated rather 
briefly as the following axiom, but we expand it afterward because 
we shall need the details later. 

Axiom 5 (Multiplication). A is a monoidal category. 

This means that it is equipped with a “multiplication” functor 0 : 
A xA^A and a “unit object” 1 that satisfy the usual associative and 
unit laws up to coherent isomorphism. Let us hrst explain “satisfying 
the laws up to isomorphism” and then discuss “coherent”. 

Associativity would mean that A0{B0C) is the same as {A0B)0C 
for any objects A, B, C (and similarly for morphisms). Associativity up 
to isomorphism means that these objects need not be equal but they 
are isomorphic, and we are given specihc isomorphisms 

ciA,B,c ■■ {A0 B) 0 C ^ A® {B 0 C) 

for all A, B,C, and furthermore these isomorphisms constitute a nat¬ 
ural transformation between functors A x A x A ^ A. 

Similarly, the requirement that the object 1 be a unit up to isomor¬ 
phism means that we are given natural isomorphisms 

\a A 0 A ^ A and pa '■ A0 1 ^ A. 

As is well-known from classical algebra, the associative law implies 
associative identities for more than three factors at a time; for example, 
if * is an associative operation, then all hve of the possible parenthe- 
sizations oia*h*c*d give the same result. The analogous result for 
categories is that any natural isomorphism a as above produces natural 
isomorphisms between any two parenthesizations of A 0 B 0 C 0 D. 
There is, however, an embarrassment of riches, as we can build, from a 
(and its inverse), several isomorphisms between such parenthesizations 
of four factors. Specihcally, the “extreme left” and “extreme right” 
parenthesizations are connected by a product of three a’s: 

{{A0B)0C)0D —)■ {A0{B0C))0D —)■ 

A0{{B0C)0 D) A0{B0{C0 D)). 

The same two parenthesizations are connected by a product of two 
other a’s: 

{{A0B)0C)0D {A0B)0{C0D) A0{B0{C0D)). 

One aspect of “coherence” is that these two transformations must 
agree, so that there is a single, well-dehned way of shifting the paren¬ 
theses from the left to the right. This requirement is often called the 
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(A0B)0 (C0B) 

((A^B)®C)®D A^{B^{C®D)) 

OlA,B,C ® Id I a® OtB,C,D 

(A 0 (B (g) C)) 0 D ---^ A (g ((B (g C) 0 D) 

Figure 1. The pentagon condition 

pentagon condition, because the diagram exhibiting these two transfor¬ 
mations together has the shape of a pentagon. In this connection, the 
hrst composition, involving three morphisms, is sometimes called the 
“long side” of the pentagon, and the second composition is the “short 
side”. In Figured! the short side is the top of the pentagon while the 
long side contains the vertical sides and the bottom. 

Another aspect of coherence is that two ways of simplifying [A g 
1) g 5 should agree, namely pa g Ib and 

(A g 1) g 5 A g (1 g 5) A^B. 

It is easy to think of other compositions of a’s, A’s, and p’s that 
should agree, for example the many ways of connecting different paren- 
thesizations of hve or more factors. Fortunately, all of these require¬ 
ments can be deduced from the two that we have exhibited here. This 
is Mac Lane’s coherence theorem, and we refer to Chapter VII of [7] 
for its precise statement, its proof, and additional information about 
monoidal categories. 

The pentagon condition will play a major role in the rest of this 
paper, because the associativity isomorphism a is often nontrivial and 
of considerable interest. The unit isomorphisms A and p, on the other 
hand, will play essentially no role, because one can safely identify IgX 
and X g 1 with X and take Xx = Px = Ix for all X. From now on, 
we will make these simplifying identihcations. 

The idea that g represents combining two anyons (or two systems 
of anyons) into a single system suggests that this operation should be 
commutative, i.e., that XgV should be naturally isomorphic to Y gX. 
The next axiom postulates the existence of such an isomorphism, with 
good behavior in connection with the associativity isomorphism a. 

Axiom 6 (Braiding) . The monoidal structure on A is equipped with 
a braiding, i.e., a natural isomorphism axx ■ XgV —)> Y gX subject to 
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two requirements, first that the following two composite isomorphisms 
be equal: 

A® {B®C) {B®C)®A'^^"' B®{C®A) 

and 

{A®B)®C {B®A)®C B®{A®C) B®{C®A), 

and, second, the analogous equality with each ax,Y replaced with 
C7y,x"^- 


A®(B® cf^^^(B ®C)®A 



(B®A)® C -—^ B®{A®C) 


Figure 2. The hexagon condition 

Recall, from Section [221 that anyons inhabit two-dimensional space 
and therefore, when two of them are interchanged, it is necessary to 
keep track of how they move around each other. A clockwise rotation 
by TT around the midpoint between them is not the same as, nor even 
deformable to, a counterclockwise rotation. So we should describe ax,Y 
not merely as switching X with Y but as doing so in a counterclockwise 
direction. The choice of direction here is a matter of convention; (7y,x~^ 
is then the clockwise rotation achieving the same interchange. Thus, 
we expect that, in general, crx,Y 7^ (If these two were always 

equal, then we would have a symmetric monoidal category rather than 
a braided one.) 

A useful picture, often used in connection with braiding, is to imagine 
the factors in a (8)-product as being lined up from left to right. Then 
the counterclockwise interchange (Tx,y amounts to moving X from the 
left of Y to the right of Y by passing X in front of Y. (Ty,x~^ also 
moves X from the left to the right of F, but it does so by passing X 
behind Y. 

The equality of the two composite morphisms in the dehnition of 
braiding is called the hexagon condition (Figure [2]). In terms of moving 
anyons around each other, it expresses the fact that moving A past 
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i? 0 C* by passing A in front of -B 0 C* is equivalent to first passing A in 
front of B and then passing A in front of C. The hexagon condition for 
crY^x~^ has a similar pictorial description with “in front of” replaced 
with “behind”. 

The last axiom in this subsection relates the multiplicative structure 
discussed here with the additive structure from the preceding subsec¬ 
tion. 

Axiom 7 (Additive-Multiplicative). (1) The monoidal unit 1 
is simple. 

(2) The product operation 0 is bilinear on morphisms. 

In more detail, item (2) here means that the function 

Hom(A, B) X Hom(C, D) — )■ Hom(A ® C^B ® D) 

given by the functoriality of 0 is bilinear with respect to the C-vector 
space structures of the Hom-sets. It follows from this, via results in [21 
Section 2.4], that 0 distributes over © on objects, i.e., that A0(F©Z) 
is canonically isomorphic to {X 0 F) © (A 0 Z). 

3.3. Duals, Twists, and Modularity. In this subsection, we collect 
some additional axioms to complete the dehnition of modular tensor 
categories. These axioms will not play a role in the computations we 
do later. We list them for the sake of completeness, but we make only 
a few comments about them and refer the reader to [9l Sections 4.3, 
4.5, and 4.7] for more thorough explanations. 

Axiom 8 (Antiparticles). For each object X of A, there is a dual 
object X*, and there are two morphisms ix ^ ^ X ^ X* and ex '■ 
X* 0 A —)■ 1, such that the compositions 

A* A* 0 A 0 A* A* 

and 

A A 0 A* 0 A X 

are equal to the identity morphisms Ix* and lx, respectively. Further¬ 
more, dualization commutes with 0 and © and preserves 1 and 0. 

For the sake of readability, we have exhibited the compositions in 
this axiom without the parentheses and associativity isomorphisms that 
technically should be there. We follow the same convention for iterated 
0 below. 

The intention behind this axiom is that, if A represents some par¬ 
ticle, then X* represents its antiparticle. The morphism ix represents 
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creation of a particle-antiparticle pair from the vacnum, and ex repre¬ 
sents annihilation of such a pair. 

The operation of dualization becomes a contravariant functor from 
A to itself if one dehnes the dual f* of a morphism f : X ^ Y to be 
the composite 

Y* r* 0 X 0 X* r* 0 r 0 X* x*. 

Axiom 9 (Rotations). There is a natural isomorphism 6 with com¬ 
ponents : X —)■ X** respecting the monoidal structure and duality 
in the sense that 

<^1 = h, = Sx ^ Sy, and 6x* = {^x*) 

By combining these 6 isomorphisms with the morphisms i and e from 
duality, one can obtain isomorphisms X —)■ X that represent twisting 
an anyon by 27r; see [9], Section 4.5] for details. 

Monoidal categories satisfying the “Antiparticles” axiom are called 
rigid, and those that also satisfy the “Rotations” axiom are called rib¬ 
bon categories. 

Axiom 10 (Modularity). For any two simple objects X and Y, let 
Sx,Y : 1 ^ 1 be the morphism 

1 = 10 1 X 0 X* 0 X 0 X* X 0 X 0 X* 0 X* 

X 0 X* 0 X 0 X* 0 X* 0 X” 0 Y* 

ex*0eY* H H H 

- > 10 1 = 1 . 

Since Hom(l, 1) = C, these morphisms sx,y constitute a matrix of 
complex numbers, with rows and columns indexed by the isomorphism 
classes of simple objects. This matrix is required to be invertible. 

Notice that, if A were not merely braided but symmetric, then the 
cr’s and the (T“^’s in this composite would cancel out, and we would 
have sx,Y = ixW where tx is the composite 

1 ^ X 0 X* X** 0 X* ^ 1, 

and similarly for ty. Thus, the matrix s described in the modularity 
axiom would be the product of a column vector by a row vector (in this 
order). Such a matrix has rank at most 1. By requiring this matrix to 
be invertible, the axiom says that, as far as the rank of this matrix is 
concerned, the braiding is as far as possible from being symmetric. 
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4. Yoneda simplification 

In this section, we point out a simplification of the additive structure 
of A, based on Yoneda’s Lemma. That lemma (see [71 Section HI.2]) 
says roughly that an object in any category is determined, up to iso¬ 
morphism, by the morphisms into it. More precisely, any category C 
is equivalent to a full subcategory of the category C of contravariant 
functors from C to the category of sets0 Under this equivalence, any 
object Y of C corresponds to the functor Hom(—,X), i.e., the functor 
sending each object U of C to the set of morphisms U ^ X and sending 
each morphism / : U —?■ U to the operation Hom(U, Y) —)■ Hom([/, Y) 
of composition with /. 

In the case of our category M, we can greatly simplify A while still 
maintaining the Yoneda equivalence. In the first place, since every 
object U of M is a finite sum, and thus in particular a coproduct, of 
simple objects, U = 0ieF S'j, morphisms U —?■ Y amount to F-indexed 
families of morphisms Sj —?■ Y. More precisely, any f : U ^ X 
is determined by the composite morphisms f o Uj : Sj X, and, 
conversely, any family of morphisms gj : Sj ^ X arises in this way 
from a unique morphism U —?■ Y. Thus, A is equivalent to a full 
subcategory of the category S of set-valued functors on the category S 
of simple objects in A. 

Up to equivalence, we need not use all the simple objects; it suf¬ 
fices to have at least one representative from each isomorphism class 
of simple objects. So we can replace the S of the preceding paragraph 
by a skeleton of it, i.e., a full subcategory iSq consisting of just one 
representative per isomorphism class. 

The structure of this new, skeletal iSq admits, thanks to the finiteness 
axiom and Proposition 0] the following description. There are finitely 
many objects. The morphisms from any object to itself form a copy of 
C. If U and V are distinct objects, then the only morphism from U to 
V is zero. 

As a result, the Yoneda embedding, simplified as above, sends each 
object Y of M to a finite family of vector spaces, indexed by the simple 
objects U in Sq, namely the vector spaces Hom(U, Y). Furthermore, 
the morphisms Y —)■ Y in M are given by arbitrary families of lin¬ 
ear maps gu '■ Hom(f/, Y) —)■ Hom(f/, Y) between corresponding vector 
spaces. The reason for “arbitrary” is that, because of the paucity of 


'^There are set-theoretic issues if C is a proper class rather than a set, but these 
issues need not concern us here. The finiteness conditions imposed on our anyon 
category A ensure that it is equivalent to a small, i.e., set-sized, category. 
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morphisms in iSq, all such families automatically satisfy the commuta¬ 
tivity conditions required in order to be natural transformations and 
thus to be morphisms in the functor category Sq. 

Summarizing, we have that, up to equivalence of categories, A can be 
described as the category whose objects (resp. morphisms) are families 
of hnite-dimensional vector spaces (resp. linear maps), indexed by the 
objects of Sq. Furthermore, it is easy to check that sums in A are 
given, via this equivalence, by direct sums of vector spaces. 

In other words, the additive structure of A is trivial. The interesting 
structure is the monoidal structure, and this can be quite complicated. 
In particular, the associativity isomorphisms a and the braiding iso¬ 
morphisms a, though given (like any morphisms) by linear maps, need 
not have a particularly simple structure. 

The analysis of the multiplicative structure of A can be facilitated 
by taking advantage of the semisimplicity of A and the fact that <8) 
distributes over ©. If we know how (8) acts on simple objects, dis- 
tributivity determines how it acts on sums of simple objects, and, by 
semisimplicity, those are all the objects. Moreover, because the asso¬ 
ciativity and braiding isomorphisms are natural, and thus in particular 
commute with the injection and projection morphisms of sums, the 
behavior of these isomorphisms on arbitrary objects is determined by 
their behavior on simple objects. Better yet, the pentagon and hexagon 
conditions will be satished in general as soon as they are satished for 
simple objects. 

Thus, the additive and multiplicative structure of A can be com¬ 
pletely described by giving 

(1) a complete list of non-isomorphic simple objects (including the 
unit or vacuum 1), 

(2) for each pair of objects in this list, their ©-product, expressed 
as a sum of objects from the list, 

(3) the associativity isomorphisms ax,Y,z for aU X, Y, Z in the list, 
and 

(4) the braiding isomorphisms ax,Y for all X, Y in the list, 

subject to the pentagon and hexagon conditions. 

We shall not be concerned here with duality and ribbon structure, 
but it could also be reduced to a consideration of the simple cases. 

Often, items (1) and (2) here determine or at least greatly constrain 
items (3) and (4) via the pentagon and hexagon conditions. One such 
situation is the subject of the next section. Other examples, both of 
strong constraints on (3) and (4) and of weak constraints can be found 
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5. Fibonacci anyons 

5.1. Definition and Additive Structure. In this section, we con¬ 
sider the special case of Fibonacci anyons. These are defined by spec¬ 
ifying the category A as follows. There are just two simple objects, 1 
(the vacuum, the unit for 0) and r. Each is its own dual. (Recall that 
Axiom [8] requires each object to have a dual; dualization is additive, so 
we need only specify the duals of the simple objects.) The monoidal 
structure is given byr(g)r = l©r (plus the fact that 1 is the unit, so 
l(8)r = r(8)l = r and 1 © 1 = 1). 

The terminology “Fibonacci anyon” comes from the fact, easily ver¬ 
ified using the distributivity of (8) over ©, that iteration of <8) gives 
T^n = • 1 © /„ • r, where the /’s are the Fibonacci numbers defined 

by the recursion /_i = 1, /g = 0, and fn+i = fn + fn-i- Here and 
below, we use the notation k ■ S to mean the sum of k copies of the 
object S of A. (The notation makes sense for arbitrary objects S, but 
we shall need it only for simple S.) 

As explained in Section 01 we can identify the category A with 
the category of pairs of finite-dimensional complex vec¬ 

tor spaces. Explicitly, an object X is identified with the pair 
(Hom(l, A), Hom(r, X)). In particular, the unit 1 in A is identified 
with (C, 0), and r is identified with (0,C). This identification respects 
the additive structure: © in A corresponds to componentwise direct 
sum of pairs of vector spaces. 

5.2. Tensor Products. The multiplicative structure of A, on the 
other hand, is quite far from componentwise tensor product of vec¬ 
tor spaces, as the latter would make r © r = r (because C © C = C). 
Our goal in the rest of this paper is to determine the multiplicative 
structure in terms of pairs of vector spaces. 

The equation r"' = /„_i-l©/„-r mentioned above already determines 
that structure as far as the objects are concerned, but there remains 
much to be said about the morphisms. 

A morphism from one pair of vector spaces (Fi, F-) to another such 
pair (VFi, 14A) is a pair of linear transformations (mi : Vi —)■ Wi,mT- : 
Vr —> Wr). We can think of it as a pair of matrices, provided we fix 
bases for all the vector spaces involved here. 

The choice of bases involves considerable arbitrariness, but there is a 
(somewhat) helpful guiding principle, namely that, if we have already 
chosen bases for two vector spaces, then the union of those bases serves 
naturally as a basis for the direct sum of those vector spaces. Some 
caution is needed, though, because the same vector space can arise as 
a direct sum in several ways and can thus have several equally natural 
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bases. Indeed, much of our work below will be finding the transforma¬ 
tions that relate such bases. 

The guiding principle tells us nothing about choosing bases for the 
one-dimensional spaces Vi and K- in the pairs 1 = (Vi,0) and r = 
(0, V-r). There isn’t even any non-zero morphism between these simple 
objects to suggest a correlation between the choice of bases. Nor do 
we get canonical bases here by evaluating compound expressions that 
fuse to r or to 1 or to a sum of these. So we might as well identify 
these one-dimensional spaces with C and use the number 1 as the basis 
vector in both of them. 

Then r(8)r = l©r = (C, C) already has a basis for each of the two 
vector spaces. Let us turn to the triple product 

r (g) (r (8) r) = r © (1 © r) = (r © 1) © (r © r) = r © (1 © r) = 1 ■ 1 © 2 • r. 

As a pair of vector spaces, it is isomorphic to (C, C^), but we have 
some additional information about it, namely that it was obtained as 
the sum of r © 1 = r and r © r = 1 © r. Our guiding principle thus 
suggests choosing a basis in that respects this sum decomposition. 
That is, one of the basis vectors in should come from the hrst r and 
the other should come from the second summand, 1 © r. 

Consider, however, the analogous computation with the other way 
of parenthesizing the triple product: 

(r © r) © r = (1 © r) © r = (1 © r) © (r © r) = r © (1 © r) = 1 • 1 © 2 • r. 

It also leads to the pair of vector spaces (C,C^), and it also provides 
a suggestion for a basis of C^. There is, however, no guarantee that 
this suggestion agrees with the one in the preceding paragraph. We 
shall see below that the two suggestions are actually guaranteed to 
disagree. We have two bases for C^, and there will be a non-trivial 
matrix transforming the one into the other. We shall hnd that this 
matrix is almost uniquely determined. 

There could, a priori, have also been two different natural bases 
for the hrst component C in although we shall see that, in this 
particular situation, they coincide. 

These basis transformation matrices, relating the bases that arise 
from r © (r © r) and from (r © r) © r, amount to the associativity 
isomorphism cIt-.t.t in the dehnition of the monoidal category A. 

Recall from Section |4] that all the associativity isomorphisms of A 
are determined by those with simple objects as subscripts. One of 
these is the Q;r,T,T mentioned just above; the others involve one or more 
I’s in the subscript. Fortunately, all those others are identity maps, 
thanks to the identihcation of 1 © X and X © 1 with X. So the entire 


24 


ANDREAS BLASS AND YURI GUREVICH 


associativity structure of A comes down to two matrices, a 2 x 2 matrix 
relating the two bases for and a number (a 1 x 1 matrix) relating 
the two bases for C. These matrices are subject to the constraint given 
by the pentagon condition (Figure [T]). Below, we shall calculate that 
constraint explicitly. It will almost uniquely determine a. 

We shall also calculate the constraint imposed by the hexagon condi¬ 
tion on the braiding isomorphisms a (Figure [2]). Again, the only com¬ 
ponent that needs to be computed is The components where at 
least one subscript is 1 are trivial, and the components with non-simple 
objects as subscripts reduce, by distributivity, to ones with simple sub¬ 
scripts. 


5.3. Notation for Basis Vectors. In order to compute the isomor¬ 
phisms a!T-,T,r and for Fibonacci anyons, we shall view them as 
matrices, using suitable bases for the relevant vector spaces, and we 
shall calculate the constraints imposed on those matrices by the pen¬ 
tagon and hexagon conditions. We begin by setting up a convenient 
notation for those bases. 

The domains and codomains of the morphisms under consideration 
are obtained from r and 1 by iterated 0. We must, of course, be care¬ 
ful about the parenthesization of such (8)-products because, as we saw 
above, different parenthesizations can lead to different bases; indeed, 
contains exactly the information about how two such bases are 
related. 

In general, given a parenthesized ^-product of r’s and I’s, we can use 
the dehning equations for Fibonacci anyons, particularly t0t = l©r, 
and the distributivity of 0 over ©, to convert the given product into a 
sum of r’s and I’s. Each summand in that sum arises from the original 
product as a result of certain choices of 1 or r when expanding some 
occurrences of r 0 r. 

For example, in the equation 


t0{t0t) = r©(l©r) = (r(8)l)©(r(8)r) = r©(l©r) = l-l©2-r 


considered above, the summand 1 at the right end of the equation arose 
from the r 0 {r 0 t) at the left end by hrst choosing the summand r 
in the evaluation of (r 0 r) at the hrst step in the equation, and then, 
after applying the distributive law at the second step, choosing the 
summand 1 in the evaluation of r (8) r at the third step. These choices 
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can be visualized as the tree 



or, in a more compressed notation, 


(r • (r • r)). 

1 T 


Here the three r’s and the parentheses describe the {8)-product r(8) (r® 
r) that we began with, and the symbols under the dots indicate the 
choice of summand at each step. The inner ■ indicates that, from the 

r 

evaluation of the inner r <8) r = 1 © r, we chose the r summand. After 
applying distributivity, that leads us to r © r, from which, as indicated 
by the outer we chose the summand 1. 

The other possible choices during the same evaluation would be writ¬ 
ten 

(r ■ (r ■ r)) and (r • (r • r)) 

T T r 1 

and depicted by the trees 



The hrst of these indicates that, as before, we chose the r summand 
when evaluating the inner ©, obtaining, when distributivity is applied, 
the summand r © r = 1 © r, but then we chose the r rather than 
the 1. The second indicates that, when evaluating the inner r © r, 
we chose the summand 1, so that, after applying distributivity, we got 
r © 1. Here, there is no choice remaining to be made; r © 1 is simply 
r. Nevertheless, we write r under the outer dot and at the root of the 
tree, to make it obvious that the hnal result here is r. 

In what follows, we shall systematically use the compressed notation, 
but the reader can easily draw the tree diagrams. Indeed, these dia¬ 
grams are just the parse trees of the compressed notations. The trees 
can also be viewed as a sort of Feynman diagrams, depicting how the 
anyons at the leaves of the tree fuse on their way to the root. 
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In our notation, we write a product of r’s or I’s with r’s or I’s 
also under the dots, to represent specific summands (1 or r) in the 
fully distributed expansion of a ^-product of r’s and I’s. To evaluate 
(X • F), hrst evaluate X and Y ; then apply ® to them; and then take 

the 1 summand in the result. To evaluate {X • F) do the same except 

r 

that you take the r summand in the result. These notations will never 
be used in situations where they would be meaningless because the 
required summand is not present in the result; that is, we never write 
{X -Y) when one of X, Y evaluates to 1 and the other to r, for then 

yields only r; and we never write {X ■ Y) when both of X, Y evaluate 

r 

to 1. As in one of the examples above, we include the subscript under 
the dot even when that subscript is forced because one of the factors 
evaluates to 1. 

Notice that our notation provides symbols, like the three examples 
above, that denote not only an object 1 or r (which can be read off 
by just looking under the outermost dot in the notation) but also a 
particular occurrence of that 1 = (C, 0) or r = (0, C) as a subspace 
(direct summand) of a specific ^-product, namely the product with 
the same factors and the same parentheses as in our notation. 

In other words, if we are given a parenthesized ^-product of I’s and 
r’s, representing the pair of vector spaces (Fi,Ft-), then by replacing 
each (g) by either ■ or •, we obtain (either a meaningless expression 

1 T 

because some required summand is absent or) a notation for a subspace 
of Vi or Vr- It denotes a subspace of Vi (resp. Vr) just in case the 
outermost (g) was replaced by • (resp. •). 

1 r 

Our notation provides names for certain summands 1 = (C, 0) or 
T = (0, C) of certain objects (Fi, F-) of the Fibonacci category A. We 
shall also use the same notation for the resulting basis vectors. That 
is, once we have a copy of, say, (C, 0) in (Fi,F-), the number 1 in C 
corresponds to some vector in Vi, and we shall use the same notation 
for this vector as for the summand. The same goes for the case of 
copies of (0, C) in (Fi, Vr); they provide vectors in F-- 

Notice that, if we begin with some parenthesized (g)-product of I’s 
and r’s, with value (Fi, F-) in A, and if we form all possible (meaning¬ 
ful) notations by replacing (g) by • or •, then the resulting vectors, as 

1 T 

described in the preceding paragraph, constitute bases for the vector 
spaces Fi and Vr- This observation is just a restatement of the fact that 
the original parenthesized g)-product is the direct sum of all the simple 
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objects obtainable by making the choices indicated by the subscripts 
in our notation. 


5.4. Associativity. Now that we have a general notation system for 
the basis vectors in parenthesized ^-products, we turn to the specific 
cases involved in associativity and the pentagon condition. 

The unique “interesting” component of associativity, Or.T.r, which we 
sometimes abbreviate as simply a, is an isomorphism from to 

r0(r{8)r), both of which are, as pairs of vector spaces, a 1-dimensional 
Vi and a 2-dimensional Vr- The first parenthesization gives a basis 
vector 

((r • r) • r) for Vi 

T 1 

and two basis vectors 

((r ■ t) ■ t) and ((r ■ r) ■ t) for Vr- 

It t t 

The second parenthesization similarly gives a basis vector 

(r • (r ■ r)) for Vi 

1 T 

and two basis vectors 


(r ■ (r • r)) and (t • (t • t)) for Vr- 

T 1 T T 

Our task is to compute the transformation a between these base^. 
This a has two components, the first relating two bases of the one¬ 
dimensional space Vi and the second relating two bases of the two- 
dimensional space Vr- These are given, respectively, by a non-zero 
number p such that 

((r ■ r) ■ r) = p(r • (r • r)) 

r 1 It 


and a non-singular matrix 


q r 
s t 


such that 


((r • r) • r) = g(r • (r ■ r)) r(r • (r • r)) 


((r . r) • r) = s(r • (r ■ r)) + t(r ■ (r • r)). 

T T T 1 T T 


®We have chosen to regard Vi and 14 as each being a single space, independent 
of the parenthesization. The different parenthesizations give (possibly) different 
bases for these spaces. An alternative view is that each parenthesization gives its 
own Vi and 14, isomorphic to C and respectively, with their standard bases, 
while a gives an isomorphism between the two Vi’s and an isomorphism between 
the two Vt’s. The two viewpoints are easily intertranslatable and the computations 
that follow would be the same in either picture. 
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Here “non-zero” for p and “non-singular” for the matrix embody the 
requirement that a is an isomorphism. 

We shall now investigate the constraints imposed on p, q, r, s, t by 
the pentagon condition. 


(r (8) r) (g) (r (g) r) 


((r <g) r) (g) r) g) r 


arrr® It 



{t ®{t 


a 


T,T^T,T 


r (g) (r (g) pr g r)) 

It ® Ur,T,T 

r ® ((r g r) ® r) 


That condition involves the g-product of four r’s, parenthesized in 
hve ways, and we shall need to consider the natural bases for all hve 
parenthesizations. Since = (C^, C^), each parenthesization will give 
two vectors as a basis for the 1 component and three as a basis for the 
r component. We begin by considering the r components, whose bases 
are displayed below. (There is no signihcance to the chosen ordering of 
the hve bases, nor the ordering of the three vectors within each basis.) 


(((r;r) ■ r) ;r) 
((r • (r;r)) ;r) 
(r • (r • (r • r))) 

r 1 r 


(((r • r) ;r) ;r) 
((r;(r ■r));r) 
(r • (r ■ (r ■ r))) 

r r 1 


(((r ■ r) • r) ■ r) 
((r;(r;r));r) 
(r ■ (r • (r • r))) 

T T T 


Each row in this picture is a basis for the 3-dimensional W; specihcally, 
it is the basis arising from the same parenthesization ofrgrgrgr 
as the parenthesization in our notation. 

When writing transformation matrices between these bases, we must 
regard each basis as given in a specihc order, because rows of a matrix 
come in an order. We (arbitrarily) choose the orders in which the bases 
are displayed above. 

The hve isomorphisms that appear in the pentagon condition amount 
to hve transformations between these bases. Let us consider these one 
at a time, beginning with the one connecting the hrst two bases in the 
table. Here we are dealing with the isomorphism 


ar(g)r,r,T : (((u g r)g)r g r) ^ ((r g r) g (r g r)). 

The hrst subscript of this a, namely r <S)T, can be decomposed as the 
sum l©r, and the naturality of a then implies that ar^r,T,T is the direct 
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sum of ai,T,r and aT,T,r- The first of these two summands is the identity, 
like all associativity isomorphisms where one of the three factors is 1. 

r'' 


The second summand is given by our matrix 


t 


As a result, we 


hnd that the transformation ar®T,T,T connecting the hrst two bases in 
our list is (taking into account the order in which the basis vectors are 
listed) 


( 0 g r\ 
1 0 0 1 
Q s t) 


In this matrix, the 1 in the (2,1) position and the two zeros in its row 
arise from the fact that the identity map ai,T-,r sends the second vector 
in our hrst basis to the hrst vector in the second basis. Had we listed 
{{{t ■ t) ■ t) ■ t) hrst rather than second in our hrst basis, the matrix for 

1 r T 

ar^T-,T,T would have been a block diagonal matrix with 1 in the upper 
left corner. 

An exactly analogous computation gives the isomorphism between 
the second and the last bases in our list: 


( g 0 r\ 
0 1 0 1 
s 0 tj 


Multipying these two matrices, we get the transformation from the 
hrst basis (parenthesized to the left) to the last (parenthesized to the 
right) that corresponds to the “short” side of the pentagon (two mor- 
phisms, across the top of the diagram). This product is 


( rs q rt\ 

q Or. 
st s j 

Turning to the long side of the pentagon (three morphisms), we hnd 
that the middle one, corresponding to rows 3 and 4 in our list of bases 
and to the bottom of the diagram, is quite analogous to the two that 
we have already computed. It is 

( g 0 r\ 

0 1 d 1 • 

s Q tj 

The remaining two isomorphisms for the long side of the pentagon 
(the vertical arrows in the diagram) are a bit diherent, as they involve 
a’s on three of the four factors and an identity map on the remaining 
factor. Let us consider aT-,T,T ® A, which connects the hrst basis in our 
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list to the third. In effect, this ignores the rightmost factor and acts 
like a on the first three factors. In other words, it is given by the same 
matrix as the transformation from the basis 


((^ • ((^: ■ t) ((^ • 

r 1 It t t 

to the basis 

(r -(r-r)) (r ■ (r ■ r)) (r • (r • r)). 

It r 1 T T 

Notice that, in each of these bases the first element is in the Vi com¬ 
ponent, so that component of a, namely p, enters the picture. Indeed, 
the matrix connecting these bases is 

/p 0 0\ 

ar,T,T (2) /r = 0 q r . 

\0 s tj 


Similarly, the remaining isomorphism on the long side of the penta¬ 
gon is also 

^p 0 0 ^ 


Ot.-t.-t = 


0 q r 
.0 s t 


Multiplying the three matrices for the long side of the pentagon, and 
equating, as the pentagon condition requires, the resulting product to 
the product that we obtained for the short side of the pentagon, we 
have 


f p'^q pfg 

prs 
, vst 


prt 

q^ + rst qr -|- 
ns rs 


V / rs 


rt\ 

= g 

0 

r 

/ \st 

s 



This is the Vr part of the pentagon condition. Before turning to the Vi 
part, let us extract as much information as possible from the matrix 
equation that we have just derived. 

Suppose, toward a contradiction, that p 7 ^ 1. Then the (1,3) and 
(3,1) components of our matrix equation give rt = st = 0, so either 
r = s = 0 ort = 0 . Ifr = s = 0 , then the ( 1 , 2 ) component of 
the matrix equation gives that g = 0 also, but this contradicts the 

^ is non-singular. There remains the case that t = 0. 

Then the (2,2) component says g = 0, the (2,3) component says r = 0, 

g r 


fact that 


and we again contradict the non-singularity of 


t 


So we have 


contradictions in all cases if p 7 ^ 1 . 

So p = 1. Now the (1,1) entry of the matrix equation gives g = rs. 
Substituting that into the (2,2) component, we get g(g + t) = 0, so 
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either g = 0 or g = —t. The hrst of these options leads, via the (1,2) 
entry, to rs = 0 and thus to a contradiction to non-singularity, as 
before. Therefore g = —t. 

From the (2,3) and (3,2) entries, we get that (g -|- t^)r = r and 
(g -|- t^)s = s. We cannot have both r = 0 and s = 0, as that would 
give g = 0 in the ( 1 , 2 ) entry and contradict non-singularity. So we 
must have g -|- = 1. In view of g = —t, this means g^ -|- g — 1 = 0 and 


therefore 


-1 ± x/5 


This evaluation of g and t, together with the earlier results 


p = 1 and rs = g. 


satisfy, as one easily checks, the entire matrix equation above. The 
least trivial item to check is the (3,3) entry, rs + t^ = which, in view 
of the equations above, becomes g — g^ = g^, i.e., 0 = g(g^ -|- g — 1 ), and 
this is true because g was obtained as a solution of g^ -|- g — 1 = 0 . 

All of the preceding calculation was based on the W component of 
we still have the Vi component of the pentagon equation to work 
out. Again, we have a list of hve bases, now for a 2-dimensional space, 
as follows. 


-^ t )-^^ t )-^ t ) (((r^r) • r) • r) 

((^;^); (^;^)) 

((r;(r ■ r))-r) ((r ; (r ; r)) • r) 

(^; ((^;^);^)) 

(r • (r • (r • r))) (r • (r • (r ■ r))) 

1 T 1 1 r T 


Computations analogous to (but shorter than) the earlier ones give, 
for the short side of the pentagon. 



So the product for the short side is simply 
we get 




and 
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Equating the product of the long side and the product of the short 
side, we get 

1 0 \ + prs qr + ptr 

0 p^j \qs + pts rs + pt^ 

This matrix equation is automatically satished because of the equa¬ 
tions that we had already derived from the component of the pen¬ 
tagon condition. So there is no new information in the Vi component. 

We can, however, get some additional information if we impose the 
requirement that the associativity isomorphisms be unitary transfor¬ 
mations. This amounts to requiring the vector spaces of morphisms 
Hom(X, Y) to be Hilbert spaces and requiring our natural bases for 
them to be orthonormal. 

Unitarity tells us nothing new about p, since we already know p = I, 

gives the equations 

+ \r\^ = = 1 and q{s — r) = q{s — f) = 0 , 

where bars denote complex conjugation and where we used the fact 
that q is real. So s = f and, since rs = q, we get hrst that q has to be 
positive, 

-l + \/5 
^ =- 2 -’ 

and second that 

r = ^/qe^^ and s = y/qe~^^ 

for some real 6. Thus, we hnally have, under the assumption of unitar¬ 
ity, 

(1 0 0 \ 

ttr.r.r = 0 q y/qC^^ 

Vo -q J 

with q = and 9 arbitrary. The presence of 9 here is an artifact of 

our choice of bases. If we modihed the hnal vector in each of our bases, 
((r • r) • r) in the domain of Ur.r.r and (u ■ (r • r)) in the codomain, 

T T T T 

by a phase factor then, with respect to the new bases, we would 
have 


but unitarity of 


q r 
s t 
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5.5. Braiding. We now turn to the task of computing the braiding a 
in the Fibonacci anyon category A. The only nontrivial component of 
the natural isomorphism a is ar^r, because components with a subscript 
1 are identity morphisms and components with non-simple subscripts 
reduce to direct sums of components with simple subscripts. 

The nontrivial component ar,T is an isomorphism from r $§ r = 1 © r 
to itself. Representing objects of A by pairs of vector spaces, we have 
that ar^T is an automorphism of (C,C), so it amounts to two non¬ 
zero scalars, a multiplying vectors in the hrst (1) component and b 
multiplying vectors in the second (r) component. These are subject to 
the hexagon identity, which equates the composites 


r' 


r' 




<TT,r(S)T 


r © r 


r 



(r © r) © r 



r © (r © r) 



(r © r) © r -—^ r © (r © r) 

^r,r,T 


as well as the analogous identity with cr“^ in place of a. 

Consider the first (1) component of this equation. In the bottom 
composition, the ar^r factors in the first and third morphisms must act 
on the r components so that the ©-product with R has a 1 component. 
So both of these are b. The a between them, acting on the 1 component, 
is an identity map, because our previous calculation gave p = 1. So 
the bottom of the hexagon is 5^. In the top, both of the a’s are again 
just 1. The a in the middle of that row is i-e-, fhe direct sum 

of and < 7 ,- ,-. The first of these two summands has no 1 component; 
the second does, and it is a. So the top of the hexagon is just a, and 
the hexagon condition reads a = b^. (The corresponding calculation 
for a~^ gives only a~^ = b~‘^, which is no new information.) 

Now consider the second (r) component of the hexagon equation. 
We do the calculation in matrix form, using the natural bases 

((r ■ t) ■ t) and ((r ■ t) ■ t) for (r © r) © r 

It t t 

and 

(r • (r • r)) and (r • (r ■ r)) for r © (r © r). 

r 1 r 1 
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With respect to these bases, is given by 

earlier. Both (Tr,r <H) It and It <E> <7t,t are given by 


\ 

1 as computed 


a 0\ _ 0\ 

0 b) 


because in each case, acts as a on the hrst basis vector (where it 
interchanges two r’s that were combined to 1) and as b on the second 
(where it interchanges two r’s that were combined to r). Finally, aT^T^T 
is the direct sum of (Tta which is 1 and acting on the r component, 
which is 6; since that direct sum decomposition matches our choice of 

bases, (Tt^t^t is given by the matrix . Multiplying the matri¬ 

ces for each of the rows, we hnd that the hexagon identity, in the r 
component, reads 

/ q'^ + brs {q + bt)r\ _ fb'^q b^r\ 
y{q + bt)s rs + bt^ J \b^s b‘^t J ’ 

Since we know, from our associativity calculation, that r and s are 
not zero, the (1,2) and (2,1) entries of this matrix equation reduce to 
q + bt — b^ = 0, or, since t = —q, 


b^ = q{l-b). 


The (1,1) and (2,2) entries give, after we remember that rs = q and 
cancel a common factor g, 

q + b = b'^ and l + bq + b‘^ = 


The last of these equations, being quadratic in b, can be solved explic¬ 
itly: 

^ ^ -q± -4 

2 

We note that, since g = is between 0 and 1, the square root in 
the formula for b is imaginary, so the two values of b are each other’s 
complex conjugates. The product of the two values for 6 is 1, so 6 is a 
complex number of absolute value 1 with real part 

The ambiguity in the choice of b is unavoidable in this situation. 
Replacing one choice by the other just replaces a by its inverse (since 
|6| = 1), and there is nothing in the algebra of A that distinguishes the 
counterclockwise motion dehning a from the clockwise motion dehning 
To put it another way, the change from one value of b to the 
other can be exactly compensated by reflecting the orientation of the 
(2-dimensional) space in which the anyons live. 
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Although we have now computed b and thus also a = 6^, we can get a 
more useful view of these numbers by manipulating the three equations 
above that relate b to q. Solving the last one for q in terms of b, and 
substituting the result, q = into the other two equations, we 

obtain from the first equation that 

63 = - - ^ \ i.e., b^-b^ + b^-b + l = 0, 

b 

which means that —b is a primitive fifth root of unity and therefore b 
is a primitive tenth root of unity. The third equation above confirms 
that by reducing to b^ = —1. 

Among the four primitive tenth roots of unity only two, e^37r2/5^ have 
negative real parts, as b does (recall that its real part is —q/2). So we 
conclude that, up to complex conjugations, 

b = e3^*/5 therefore a = 

This completes the calculation of the braiding a for Fibonacci anyons. 

Remark 2. The multiplicative structure for Fibonacci anyons, summa¬ 
rized by the fusion rule = l0r, is perhaps the simplest nontrivial 
fusion rule. Other fusion rules have been analyzed, either by hand as 
we have done here or with computer support. The appendix of [2] sum¬ 
marizes much of what is known about specific examples. There does 
not, however, seem to be any general theory for arbitrary fusion rules. 

5.6. Fibonnaci Anyons and Qnantum Computation. In Sec¬ 
tion |2l we mentioned the hope that, by using anyons to encode qubits, 
one could use braiding to transform anyon states in various ways, 
thereby enabling quantum computation. Two anyons are not sufficient 
for this purpose, because the braid group on two strands is abelian, 
whereas quantum computation needs non-commuting unitary trans¬ 
formations. In the case of Fibonacci anyons, the computation in the 
preceding subsection shows that the braiding transformation (Jr^r is 
diagonal in a suitable basis, so it splits into one-dimensional represen¬ 
tations; this again shows its inadequacy for quantum computation. 

With three Fibonacci anyons, the situation improves dramatically. 
In a suitable basis, the transformation that braids the first two of the 
three anyons, ar^r^Iri is still diagonal. The same goes for the transfor¬ 
mation that braids the second and third anyons, but the suitable bases 
in these two cases are not the same. They differ by an associativity 
isomorphism a. More precisely, one is the conjugate of the other by 
ar,T,T- They do not commute. 

In fact, such braiding transformations suffice to approximate arbi¬ 
trary unitary transformations of the two-dimensional Hilbert space W 
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for Furthermore, using six Fibonacci anyons to code two qubits, 
one can approximate, by braiding, the so-called “controlled not” gate, 
which, in combination with one-qubit gates, is sufficient to produce 
all unitary gates for an arbitrary number of qubits; that is, it is suf- 
hcient for quantum computation. We refer to [9l Section 6] for these 
combinations of Fibonacci braidings. 
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